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ABSTRACT 

We compute quasar-galaxy and AGN-galaxy cross- correlation functions for 
samples taken from the [Veron-Cetty fc Veron ( 1998 )| catalog of quasars and 



active galaxies, using tracer galaxies taken from the Edinburgh/Durham South- 
ern Catalog. The sample of active galaxy targets shows positive correlation at 
projected separations r p < Qh^ 1 Mpc consistent with the usual power-law. On 
the other hand, we do not find a statistically significant positive quasar-galaxy 
correlation signal except in the range 3h _1 Mpc < r p < dh^Mpc where we 
find similar AGN-galaxy and quasar-galaxy correlation amplitudes. At separa- 
tions r p < 3 h^ 1 Mpc a strong decline of quasar-galaxy correlations is observed, 
suggesting a significant local influence of quasars in galaxy formation. In an at- 
tempt to reproduce the observed cross-correlation between quasars and galaxies, 
we have performed CDM cosmological hydrodynamical simulations and tested 
the viability of a scenario based on the model developed by [Silk fc Rees (1998 



In this scheme a fraction of the energy released by quasars is considered to be 
transferred into the baryonic component of the intergalactic medium in the form 
of winds. The results of the simulations suggest that the shape of the observed 
quasar-galaxy cross-correlation function could be understood in a scenario where 
a substantial amount of energy is transferred to the medium at the redshift of 
maximum quasar activity. 



Subject headings: cosmology: large-scale structure of universe — quasars: clustering- 
quasars: statistics — galaxies: clustering 
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1. Introduction 



Since the pioneer studies on the cosmological nature of quasar redshifts (eg. |Bahcall et 
al. (1969)1 , |Gunn (1971)| ) several works have provided useful insights on the environment of 
quasars such as the recent works of Smith et al.(1995) and Fischer et al. (1996) . At low 
redshifts, z < 0.3, recent results suggest that quasars are placed in similar environments 
to galaxies QSmith et al.( 1995)1) , or i n groups of 10-20 galaxies ([Fischer et al. ( 1996 )| ) . In 
this respect at high redshifts, z > 0.6, the evidence of different environments for radio-loud 
and radio-quiet quasars is compelling (|Yee fc Green (1987)| , [Boyle et al.( 1988)1 , |Yee (19901 ). 
Radio-loud quasars are mainly associated with rich Abell clusters and radio-quiet quasars 
reside either in groups or in the outskirts of clusters. It should be noted that due to the 
lack of suitable low redshift quasar samples the majority of observational works concern 
intermediate redshifts, 0.3 < z < 0.6. 

On the theoretical side the discussion of the influence of quasar activity on the formation 
and evolution of nearby galaxies has been present during the last decades. The surrounding 
galaxy environment would be responsible for triggering and fueling the quasars as first sug- 
gested by |Toomre fc Toomre (1972"]] , where interacting galaxies and mergers are efficient at 
transporting gas into the inner regions of the host galaxy. The evolution of the rate of galaxy 
interactions in clusters would allow an understanding of the observed clustering of galaxies 
around radio-loud quasars ( jStocke fc Perrenod (19817)1 ) ■ On the other hand, the large range 



of influence of a quasar may have implications for structure formation. According to |Rees 
1988)| and [Babul fc White ( 1991 )| quasars may ionize the surrounding medium up to several 



megaparsecs away. A different approach is adopted by other authors (see for instance Vbit 



|(1996)| , |Natarajan et al. (1997)| , |Silk fc Rees (1998)] ) , where a fraction of the energy released 
by quasars is transferred via energetic outflows to the gaseous medium. 

In order to shed some light on this subject we have computed the angular and pro- 
jected cross-correlation function of target quasars and AGNs from the |Veron-Cetty fc Veron 



|(1998)| Catalog, and galaxy tracers from the Edinburgh-Durham Galaxy Catalog ([Heydon- 



Pumbleton et al.(1989)| ). We have also explored the possible astrophysical effects of quasar 
activity by means of hydrodynamical cosmological numerical simulations which take into ac- 
count momentum transfer to the gas at the time of maximum quasar activity. As an attempt 
to model these effects we have considered here a simple scheme following the model of |Silk fc 
Rees (1998][. Quasar-galaxy association in the models are quantified by computing the spa- 



tial quasar-galaxy cross-correlation function and we derive the projected cross- correlations 
which may be compared to observations. 
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2. Data 

Our target sample consists of quasars and AGNs selected from the |Veron-Cetty fc Veron| 



(1998)| Catalog, hereafter VCV, which comprises angular positions, redshifts, and photomet- 
ric data for 11,358 quasars and 3,334 AGNs and provides the most recent compilation of 
quasars and AGNs. 



The Edinburgh-Durham Southern Galaxy Catalog, hereafter COSMOS Survey flHeydon- 



Pumbleton et al.(1989)| ), was used as a suitable tracer of the galaxy distribution. This survey 
provides angular positions and photographic magnitudes in the bj band for over two million 
galaxies. We have restricted our study to the region 5 < —5° given the lower quality in the 
photographic material in the northern hemisphere. We adopt a limiting apparent magnitude 
bj = 20.5 as a compromise between plate quality and depth. 

We have estimated the galaxy redshift distribution function for different apparent limit- 
ing magnitudes as derived by |Baugh fc Efstathiou ( 1993 )| (see figure 1). Taking into account 
these results we have adopted a target redshift cutoff z = 0.35 since the fraction of COS- 
MOS galaxies beyond this redshift is less than 2% for limiting magnitude bj = 20.5. We 
have excluded from our analysis targets with z < 0.05. This avoids a few nearby objects 
for which fixed projected spatial bins correspond to very large angular areas dominating 
projected correlations. With the above restrictions the resulting sample of targets consists 
in 103 quasars and 226 AGNs in the region a < 4M0™ and a > 21 1 ?, and -60° < 5 < -5°. 
The tracer sample comprises 1,982,847 galaxies brighter than bj = 20.5 in the same area. 



3. Analysis 

Our analysis consists of the computation of the mean excess with respect to the average 
number density of COSMOS galaxies as a function of angular separation w(9)), and projected 
distance (W(r p )) from VCV targets. 



w{9) 



(N r (9)) 
(Nt(r p )) 



(N t (8)), (N r (9)) ,(N t (r p )) and (N r (r p )) are the mean number of target-galaxy pairs 
(VCV-COSMOS) and random-galaxy pairs (random-COSMOS) at angular separation 9 and 
projected separation r p respectively. 
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We have considered cross-correlations for 9 > 10 arcsec and r p > 0.1 ft, Mpc to 
exclude the cross-correlation with the host galaxy. In order to directly compare the results 
of quasars and AGNs we have restricted the redshift distribution of AGNs to fit that of 
quasars. This new restricted sample of AGNs has 116 objects and was obtained by imposing 
to the original AGNs sample a radial gradient through a Monte-Carlo rejection algorithm that 
provides a similar redshift distribution to quasars. Therefore, both angular and projected 
cross-correlation functions can be directly compared. In figure 1 we show the corresponding 
redshift distributions of galaxies with bj < 20.5 ( [Baugh fc Efstathiou (1993]D , quasars, AGNs, 



and restricted AGNs. As can be seen, in this figure a negligible fraction of COSMOS galaxies 
are expected beyond our limiting redshift z < 0.35. The different redshift distributions of 
quasars and AGNs and the similar redshift distributions of quasars and restricted AGNs can 
also be appreciated. 

In figures 2 to 5 are shown the resulting angular and projected correlation for quasars 
(w Q (9), W Q (r p )) and the restricted AGN sample (w A (0), W A {r p )), both with COSMOS trac- 
ers at limiting magnitude bj < 20.5. In these figures error bars correspond to uncertainties 
derived from the bootstrap resampling technique developed by |Barrow et al. (1984)| with 



20 bootstrap target samples. Figure 2 shows a lack of strong quasar-galaxy angular corre- 
lations wq(9), consistent with the results derived by |Smith et al.(1995")| . On the other hand 
the AGN-tracer angular cross-correlation function wa{9) as shown in figure 3 is consistent 
with a power-law with significant correlation signal for 9 < 1,200 arcsec. 

As can be seen in figure 4 the quasar-tracer projected cross-correlation function Wq(r p ) 
shows a statistically significant positive signal for quasar-tracer projected separations in the 
range ?>h~ x Mpc < r p < Qh~ l Mpc. On smaller scales, r p < 3 h^ 1 Mpc, there is a marked 
decrease of the correlation amplitude which makes Wq(r p ) ~ 0. This behavior contrasts 
with the smooth power- law trend of Wa{t p ) shown in figure 5 for the restricted sample of 
AGN targets. Thus, a significant difference in the environment of quasars and active galaxies 
is implied by our analysis. While the latter show a correlation function similar to that of 
the general galaxy population, the projected quasar-galaxy correlation function, Wq(r p ), 
indicates a significant departure from a power law at small scales. 

By comparing the projected correlation functions Wq{t p ) and WA{r p ) it can be seen their 
similar correlation amplitude on scales 3h~ l Mpc < r p < Qh^Mpc. This result indicates that 
quasars and AGNs reside in similar global galaxy density enhancements and suggesting that 
the effective influence of quasars on galaxy formation is restricted to scales < 3h~ 1 Mpc. 
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4. Numerical models 

In an attempt to analyze the influence of quasars on environment we have performed 
numerical simulations using the Hydra code, developed and kindly made available by H. 
Couchman. Initial positions and velocities of particles were generated using the Zeldovich 
approximation corresponding to the CDM power spectrum. We have adopted two different 
values of the density parameter Q = 0.2 and Q = 0.4 both with baryonic density parameter 
Qb = 0.1. Simulations corresponding to Q — 0.4 model were normalized according to Cosmic 
background temperature fluctuations as determined by measurements of the COBE satellite, 
which implies present rms mass fluctuation in a 8 h~ x Mpc sphere radius 0g = 0.75 ( |G6rski 



et al (1995)|). We have taken the same Ug value for the Q = 0.2 model in order to make them 



of comparable present mass density fluctuations and to emphasize the dependence of the 
results on the ratio of baryonic-to-total density. The computational volume corresponds to 
a periodic cube of side length 20h~ 1 Mpc measured at z = 0. We have followed the evolution 
of N = 2 x 32 3 particles for the two CDM models adopted, Q = 0.4 and Q = 0.2 with the 
same value of the Hubble constant (Hq = 65 km s _1 Mpc _1 ). The mass of the gas and the 
dark matter particles are 2.9 x 10 9 Mq and 8.59 x 10 9 Mq respectively for the Q = 0.4 model 
and 2.9 x 1O 9 M and 2.9 x 1O 9 M for the Q = 0.2 model. 

We use an analytic fit to the CDM power spectrum of the form 



2 



v(k) K k _ An(l + 2.34g) 

n ) [l + 3.89g+(16.1g) 2 + (5.46g) 3 + (6.71g) 4 ] 1/2 V 2 .34/g J K ' 



where q = (k/T)h Mpc ~\ T = Q h exp (-(0 B +y (/i/0.5)O B /O )) and Q B = 0.0125 h' 2 
is the value of the baryon density parameter imposed by nucleosynthesis theory ( Bardeen et 



al. (1986)1 an d [Sugiyama (1995)1 ). The initial conditions correspond to redshift z = 40 and 



the evolution was followed using a variable time-step resulting typical runs of 1000 steps. 

We have adopted the start of quasar activity in the models at z = 3 which coincides ap- 
proximately with the redshift of observed maximum abundance of quasars. At this redshift, 
we identify local mass density maxima and associate the quasar with the highest density 
particle. 

In order to analyze the effects of quasar energetics on galaxy formation we have applied 
a simple scheme following the model developed by Silk fc Kees (199~%) in hydrodynamical 



numerical simulations. In this plain scheme it is assumed that a fraction (~ 10%) of the 
total energy input by a quasar (~ 10 62 erg) is in the form of momentum transfer to the 
surrounding gas. This is accomplished by increasing the radial component of the velocities 
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of gas particles within a given radius at z = 3 where we assume one quasar becomes active 
in each simulation. According to |Silk fc Rees (1998)| (see also |Voit (1996)| ) this increase 



of radial velocities amounts to ~ 500A;ms 1 for particles within a sphere of radius R = 
4.6(£ 62 /fi .05) 1/5 (l + z)-^ 5 Mpc. Here, E 62 refers to the total energy in units of 10 erg 
emitted by the quasar in its lifetime (~ 10 7 yr), and f2 .05 is the baryonic density parameter 
in units of 0.05. The quasar position is chosen to coincide with the highest gas density particle 
in the simulations. It should be recalled that the simulated box is not typical regarding the 
number density of quasars since to match observed quasar abundance estimates only 1 out 
of ~ 100 such boxes should have a quasar. Our rather small box size adopted is mainly to 
improve the numerical resolution of the simulations. 

Using a standard friends-of-friends algorithm we define gaseous clumps that should be 
associated to galactic halos in the simulations. Each given initial condition is evolved with 
and without quasar activity. We take advantage of the identical initial conditions for each 
pair of simulations for a suitable analysis of the environmental effects of quasars. We compute 
the spatial quasar-galaxy correlation function in the set of simulations with quasar energy 
release and we compare the quasar-galaxy cross-correlation function with the galaxy-galaxy 
correlation functions of the set of simulations without quasars. For the latter, the centers of 
the galaxy-galaxy correlations correspond to the maximum gas density particles that defined 
the quasars at z = 3 in the former. 

The resulting spatial quasar-galaxy correlation function £q(t") corresponding to 20 pairs 
of simulations are shown in figures 6 and 7 (Q = 0.2 and Q = 0.4 respectively) which display 
the strong suppression of galaxy formation in the neighborhood of the quasars (r < 2h~ 1 Mpc) 
consistent with no correlation at these scales, and a maximum in £q(v) at r ~ 3h~ 1 Mpc 
since the power-law behavior at larger scales is not affected by the model. For comparison, 
in these figures is also shown the spatial galaxy-galaxy correlation function with center 
galaxies defined in the previous paragraph. In both figures, errorbars correspond to Poisson 
estimates, 5^q = y/N /N u where N is the number of galaxies in each bin and N u is the 
expected number in a uniform distribution. 

Comparing the spatial correlation functions of the models a stronger galaxy suppression 
at small scales is seen in the lower density (Q = 0.2) model. This is due to the higher relative 
fraction of baryons which are affected by the momentum transfer in this case. In the Q = 0.4 
model the dark matter is dominant by a factor 3 allowing the retention of more gas in the 
quasar neighborhood than in the low density models. 

In figures 8 and 9, we show the projected quasar-galaxy correlation function of the 
models W(r p ). In the small angle approximation Wq{t p ) is related to the spatial correlation 
function £g(r) by: 
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/•(JXJ 

W Q (r p ) = B ^((rj + vr 2 ) 1 / 2 )^ (4) 



with B = J2iP(yi)/(J2i(^-/yi) Io° p(x)x 2 dx) is a normalization constant which includes 
the probability p(x) that a tracer galaxy is at distance x, and yi is the distance to the i th 
quasar. We have used a Schechter function fit flSchechter (1976)]) to the luminosity function 



of COSMOS tracer galaxies, with parameters M* = —19.7, a = —1.0, a flat cosmology, and 
a linear A'-correction term K = 3z. 

The model projected quasar-galaxy correlation function on scales r p < 3 h^ 1 Mpc has 
a similar behavior than the observed cross-correlation. By inversion of equation 4 using the 
observed Wq{t p ) we derive negative values of the spatial cross-correlation function £q(V) on 
these scales, which stress the relevance of this observational effect. 

The simplicity of our adopted scheme of energy release of quasars, the small size of 
the simulated box and the restricted dynamical range impose further limitations on the 
models. Considering these shortcomings, we conclude that the results of the models analyzed 
reasonably reproduce the observations. 



5. Conclusions 

The quasar-galaxy projected cross-correlation function Wq{r v ) shows a maximum at 
~ 3h~ l Mpc and is consistent with no correlation for r p < 3h~ 1 Mpc. The inversion of such 
a projected correlation function requires a negative spatial correlation on scales < 2h~ 1 Mpc 
indicating strong effects of quasars in their local environment. The corresponding angular 
cross-correlations are roughly consistent with no signal detection. In contrast to the behavior 
of quasar-tracer correlations, we find that both angular and projected AGN-galaxy cross- 
correlation functions are consistent with the usual power-law shape. The agreement of 
Wq{t p ) and Wa(t p ) at scales 3 h~ x Mpc < r p < 6 hr x Mpc, suggests that both quasars and 
AGNs reside in similar global galaxy enhancements. 

The analysis of cosmological hydrodynamical simulations allows us to study the effects 
of energy input by the quasars to the surrounding gaseous medium. The model developed 
by [Silk fc Kees (1998)| considers a typical energy input of ~ 10 61 erg to the gas which 



develops a shell of present radius ~ AMpc originally at a speed of ~ 500 km s . The 
implementation of this scheme in our numerical simulations shows that a relative lack of 
galaxies in the neighborhood of quasars may arise naturally in these models providing support 
to the viability of this scenario. 
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Fig. 1. — Redshift distributions. The thick solid line corresponds to tracer COSMOS galaxies 
with bj < 20.5 according to |Baugh fc Etstathiou (1993)| , the thin solid line corresponds to 
quasars, dashed lines to AGNs, and dotted lines to restricted AGNs. 



Fig. 2. — Quasar-tracer angular cross-correlation function wq{9). Tracer galaxies correspond 
to COSMOS survey with bj < 20.5. 

Fig. 3. — AGN-tracer angular cross-correlation function wa{0). Tracer galaxies correspond 
to COSMOS survey with bj < 20.5. 

Fig. 4. — Quasar-tracer projected cross-correlation function Wq(t p ). Tracer galaxies corre- 
spond to COSMOS survey with bj < 20.5. 

Fig. 5. — AGN-tracer projected cross-correlation function Wa{t p ). Tracer galaxies corre- 
spond to COSMOS survey with bj < 20.5. 

Fig. 6. — Circles: spatial quasar-galaxy cross-correlation function £,q(t) in the Q = 0.2 
models. Squares correspond to the galaxy-galaxy correlation function in models with no 
quasar activity as described in the text. 

Fig. 7. — Spatial quasar-galaxy cross-correlation function ^g(r) in the Q = 0.4 numerical 
simulations. All symbols are as in figure 6. 

Fig. 8. — Projected quasar-galaxy cross-correlation function Wq{r p ) for the Q = 0.2 nu- 
merical simulations. Solid and dashed lines correspond to models with and without quasar 
activity respectively. 



Fig. 9. — Projected quasar-galaxy cross-correlation function Wq(r p ) for the Q = 0.4 numer- 
ical simulations. Lines represent the same as in figure 8. 
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